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ABSTRACT 


This  report  demonstrates,  from  experimental  studies,  that  the  classic 
critical  loads  computed  from  small  displacement  theory  are  correct  for 
perfect  shells  under  ideal  loading  conditions.  The  method  used  is  an 
implicit  rather  than  an  explicit  one.  It  makes  use  of  the  fact  that  for  a 
realistic  structure,  the  elastic  deformation  can  be  associated  with  the 
initial  displacement  from  ideal  form,  the  load  which  produces  the  motion, 
and  the  classic  load,  for  the  structure  by  the  hyperbolic  expression 


By  choosing  the  variables  to  be  6/p  and  6,  the  relationship  can  be 
presented  in  the  form  of  a  straight  line  whose  slope  corresponds  to  the 
critical  load  for  the  ideal  case.  Ihis  is,  in  essence,  the  method 
developed  by  Aiyton  and  Perry,  in  1889,  to  analyze  column  data.  The 
generalization  to  include  other  structures  was  foreseen  by  Southwell  in  his 
classic  paper  of  1932,  although  he  offered  no  proof.  In  effect,  a  general 
proof  exists  in  the  theory  of  elastic  stability  as  presented  by  Wester- 
gaard  in  1922.  However,  the  practical  application  to  shells  has,  until  now, 
not  been  made. 

By  using  this  technique,  the  behaviour  of  cylinders  under  the  action  of 
external  pressure,  torsion,  and  axial  load  and  of  the  instability  of 
spheres  and  spherical  caps  under  the  action  of  external  pressure  is 
examined.  In  the  main,  the  analyses  are  conducted  on  experimental  data 
already  published.  However,  there  are  notable  exceptions,  the  cylindrical 
shell  under  axial  load  being  a  case  in  point. 
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INTRODUCTION 


Shell  structures  are  very  common.  For  over  a  century,  engineers  have 
investigated  experimentally  and  discussed  theoretically  their  behavior 
under  various  loading  conditions.  The  extent  of  this  endeavor  is  clear 
from  the  fact  that  some  23391>2books  and  papers  were  written  in  the  field 
prior  to  1956.  Nevertheless,  many  facets  of  the  subject  are  still  unre¬ 
solved.  Indeed,  it  has  become  traditional  to  expect  wide  divergence 
between  experimental  results  and  theoretical  predictions  when  certain  prob¬ 
lems  of  stability  are  considered.  This  is  particularly  so  for  circular 
cylindrical  shells  under  axial  compression  and  spherical  shells  under 
external  pressure.  The  exact  causes  of  the  discrepancy  are  hard  to  iden¬ 
tify.  However,  one  factor  of  importance  is  easily  recognized.  Theories 
are  normally  developed  for  ideal  bodies  with  perfect  loading  conditions 
and  boundary  restraint;  whereas  tests  are  always  made  with  realistic 
structures  imperfectly  loaded  and  constrained.  Thus,  if  agreement  between 
theory  and  experiment  is  to  be  achieved,  it  will  be  through  an  implicit 
rather  than  an  explicit  process.  The  Southwell3  technique,  which  has 
already  been  demonstrated  to  have  general  applicability  in  questions  of 
elastic  stability  of  plates  and  columns,  is  such  a  process.  The  present 
paper  establishes  the  validity  of  this  technique  for  shell  problems. 

The  analytical  basis  of  the  work  is  older  than  the  Southwell  concept  itself, 
being  found  in  the  study  made  by  Westergaard5in  1922.  In  this  very  general 
analysis,  the  relationship  among  elastic  deformation,  the  initial  deviation 
from  ideal  form,  the  load  level,  the  critical  load  and  its  harmonics  was 
derived.  The  starting  point  was  the  Lagrangian  expressed  in  generalized 
coordinates.  Westergaard  considered  not  only  initial  geometric  variances 
but  also  nonhomogeneity  and  eccentricities  of  loading.  He  applied  the 
theory  of  minimum  potential  and  obtained  a  general  relationship  among  the 
actual  load  (P),  the  critical  load  for  a  particular  mode  (nPCr)>  the  ampli¬ 
tude  of  that  component  of  the  initial  irregularity  which  corresponds  to  the 
nthmode  ( n^o) >  a^d  the  appropriate  associated  component  of  elastic  deflec¬ 
tion  ( 6) . 


This  relationship  is 


6  = 


P 


P  -P 


n  cr 


6 

n  o 


(1) 


Clearly,  there  are  n  terms 
of  the  structure  under  the 


of  this  type  and  the  total  elastic  deflection 
given  load  system  is 


P  -P 
n  cr 


n6o 


(2) 


which  is,  of  course,  identical  to  the  Southwell  formulation  for  the  column. 

In  many  problems,  the  first  mode  predominates,  and,  thus,  a  single-term 
expression  of  the  type  given  in  equation  (2)  adequately  describes  the 
behavior.  Hence,  the  relationship  between  6  and  P  is  hyperbolic.  However, 
if  the  variables  are  taken  as  6/P  and  6,  the  relationship  is  linear. 
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EXPERIMENTAL  EVIDENCE 


To  illustrate  the  process,  we  chose  the  simplest  curved  structure  with  the 
most  elementary  loading;  viz,  a  circular  grch  compressed  by  a  point  load  at 
its  vertex.  Langhaar,  Boresi,  and  Carver  have  examined  this  problem  both 
theoretically  and  experimentally.  The  arch  which  they  used  was  of  10- inches 
radius  and  was  made  from  a  strip  of  aluminum  2.5  in.  wide  and  .031  in.  thick. 
It  was  simply  supported  at  its  ends.  The  load  displacement  relationship 
obtained  in  the  test  is  given  in  Figure  1.  When  this  is  plotted  in  the 
Southwell  fashion,  the  characteristic  straight  line  (Figure  l)  is  obtained. 
The  critical  load  determined  from  this  line  is  U.12  lbs,  which  is  in  excel¬ 
lent  agreement  with  the  theoretically  predicted  value  of  4.17  lbs. 


A  somewhat  more  complex  problem  is  the  behavior  of  a  tube  under  external 
pressure  (p).  This  problem  has  been  treated  analytically  by  Sturm  ,  and 


the  validity  of  his  solution  has  been  well  established  by  his  own  tests 
and  those  of  Cleaver.®  However,  from  the  viewpoint  of  this  paper,  the  most 
important  factor  in  Sturm's  analysis  lies  in  his  demonstration  that  for  a 
slightly  imperfect  cylinder,  the  initial  radial  deformation  amplitude  is 
related  to  the  subsequent  amplitude  by  the  equation 


(3) 


cr 

This  equation  is,  of  course,  identical  to  the  general  equation  derived  by 
Westergaard  and  clearly  shows  that  a  Southwell  representation  of  test  data 
should  be  applicable.  When  the  load  displacement  data  given  in  Figure  l6 
and  18  of  Sturm's  paper  are  replotted  in  the  Southwell  form,  the  straight 
lines  given  in  Figures  2  and  3  of  this  report  are  obtained.  The  critical 
loads  obtained  from  the  several  slopes  are  in  excellent  agreement  with  the 
theoretical  predictions. 


Subsequent  to  the  work  of  Sturm,  a  number'  of  experimental  studies  were  made 
on  the  stability  of  ring- reinforced  circular  cylinders  under  external 
pressure  loading,  (in  one  particular  series  of  tests,  the  shells  failed  by 
general  instability.)  It  is  with  reference  to  the  interpretation  of  these 
tests  that  the  first  practical  application„of  the  Southwell  method  to  shell 
bodies  is  recorded.  Galletly  and  Reynolds^  made  the  analysis.  In  their 
paper,  they  noted  the  applicability  of  the  Southwell  process  to  many  prob¬ 
lems  and  suggested  a  modification  for  the  case  in  question.  This  was  the 
use  of  strains  instead  of  displacements,  a  step  which  is  readily  seen 
admissable  and  which  is  analytically  demonstrated  in  Westergaard 's  paper. 

The  results  which  they  obtained  are  most  important.  For  the  first  test 
which  they  report,  observations  were  made  au  21  different  points  on  the 
body.  The  critical  loads  derived  from  these  measurements  are  in  excellent 
agreeraent  with  one  another.  They  are  listed  in  Table  I. 

It  is  seen  that  the  critical  loads  range  from  l68lb  to  178  lb,  with  a  mean 
of  172.2  lb.  These  values  compare  vslj  with  the  theoretical  load  of  173  lb 
computed  in  accordance  with  Kendrick's  theory. 

The  second  set  of  data  presented  'ry  these  researchers  is  a  comparison  of  the 


FIGURE  2.  Circular  Cylindrical  Shell  under  External 
Pressure  (Data  from  Fig  l6,  Ref  7). 
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6/P  -  M/LB 


A/W  -  IN.*/LB  *10 


FIGURE  3-  Circular  Cylindrical  Shell  Under 
External  Pressure  (Data  From 
Figure  18,  Ref  7) . 
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results  of  tests  made  on  five  similarly  constructed  shells,  whose  theoretical 
collapse  pressures  are  150  nsi  and  1000  psi.  The  ratios  of  the  critical  load 
levels  computed  from  Kendrick's  theory  and  the  values  deduced  from  the  tests 
are  given  in  Table  II.  Again,  the  comparisons  are  seen  to  be  very  good  indeed. 

11 

Donnell  in  his  classic  paper  on  the  Southwell  method  appears  to  have  been 
the  first  to  consider  whether  or  not  the  process  could  be  applied  to  the 
thin-walled  circular  cylinder  under  axial  load..  He  examined  the  problem 
analytically,  taking  as  his  starting  point  his  approximate  solution  for  the 
instability  of  a  cylinder  under  axial  compression  in  which  he  had  considered 
initial  deviations  from  cylindrical  shape. 


He  chose  the  deflection  function  to  be  given  by 

w  =  W  sin  5—  sin  ——  +  W(W+2w' )  ~ 

L  L  Sir 


cos 


2TTx 


(«0 


The  reason  for  this  choice  lay  in  the  fact  that  ;he  second  tenn  was  most 
useful  in  annulling  some  of  the  second-order  extension  with  first-order 
extensions  and,  thus,  was  effective  in  reducing  the  internal  energy.  This 
initial  deviation  from  cylindrical  form  was  taken  to  be  given  by 


w'  = 


W 

w" 


w 


(5) 


From  this  basis,  he  computed  that  the  relationship  among  the  load,  the 
critical  load,  and  the  initial  and  subsequent  deformations  would  be  given 
by 


1  + 


^  W  P 

p  =  vr7- v*  cr 


l  + 


ttV 

8l~ 


( W+2W ' )  (W+W') 


(W+2W)2  (2W*W) 
WrW 


(6) 


where  P  is  the  critical  load  for  the  perfect  shell  and  is  given  by  the 
normal  Equation 


i  .  nut 


(7) 


It  is  clear  that  if  Donnell's  large  displacement  derivation  is  true,  there 
must  be  some  classes  of  shells  for  which  the  multiplying  factor  on  the 
Southwell  form  is  not  unity.  But  it  is  apparent  that  if  W  and  W'  are  smell 
compared  with  2L  /TrR  or  when  W  is  very  small  compared  with  W ,  the  South- 
well  plot  should  yield  accurate  results. 


In  view  of  Donnell's  analytical  development  and.  the  acute  interest  in  circu¬ 
lar  cylinders  under  axial  compiession,  it  might  be  thought  that  the  question 
would  have  been  experimentally  resolved.  Unfortunately,  at  the  time  that 
Donnell  was  studying  this  problem,  it  was  not  technically  feasible  to  use 
noncontacting  displacement  transducers  to  measure  the  radial  movement  of  the 
wall  of  an  axially  compressed  shell  model  with  the  degree  of  accuracy 


6 


7 


required  to  verify  the  analysis.  Thus,  Donnell  .ms  constrained  to 
theoretical  observation  only. 

i  12 

Flugge,  in  his  text  on  "Stresses  in  Shells^  appears  to  be  the  next  author 
to  refer  application  of  the  Southwell  plot  to  test  results  on  shell  bodies. 
In  his  discussion,  he  emphasizes  that  more  than  one  mode  may  be  involved  in 
the  total  deformation;  thus,  since  in  many  shells  there  exist  a  number  of 
critical  loads  immediately  above  the  lowest  one,  several  terms  of  the  dis¬ 
placement  series  may  grow  to  infinity  together.  In  this  detail,  then,  the 
cylinder  under  axial  compression  differs  from  the  column  for  which  the 
eigenvalues  are  significantly  separated. 

However,  he  states  that  if  care  is  taken  to  measure  a  displacement  which  is 
large  for  the  lowest  buckling  load  but  small  for  all  others,  then  the 
Southwell  technique  is  applicable.  He  substantiated  his  remark  by  refer¬ 
ence  to  unpublished  experiments  conducted  by  Kroram  end  Flugge  in  the  early 
19*4-0' s.  These  tests  were  made  on  cylinders  under  various  combinations  of 
axial  compression  and  torsion.  Professor  Flttgge  has  remarked  that  Southwell 
plots  were  obtained  only  when  the  predominant  loading  was  torsional. 

In  our  research,  we  have  found  that  if  well-made  thin-walled  shells  are 
tested  with  extreme  care  to  ensure  uniformity  of  load  distribution,  and  if 
noncontacting  probes  of  high  resolution  are  used  to  determine  the  wall 
motion,  load-displacement  relationships  of  the  hyperbolic  type  are  obtained. 
Data  for  typical  tests  are  given  in  Figures  k  and  5.  Both  these  observa¬ 
tions.,  were  made  on  the  same  shell.  The  data  have  been  plotted  in  the  Lund- 
quist  ^  manner.  The  agreement  between  the  classic  critical  loads  and  the 
values  computed  for  these  plots  is  good  (maximum  error,  1*98$).  The 
cylinder  used  was  made  from  a  thick-walled  tube  of  aluminum  by  machining 
and  was  circular  to  within  l/l,000  in.  and  was  uniform  in  thickness  to 
better  than  l/l0,000  in.  It  had  an  R/t  of  335  and  an^/D  of  I.87.  The 
test  was  made  in  which  loading  was  via  a  thermal  ram.  Deflection 
measurements  were  made  using  a^Fotonic  sensor,  an  electro-optical  measur¬ 
ing  device  of  high  resolution.  The  test  is  fully  described  in  reference 
l6. 

In  the  case  of  torsional  loading,  it  is  much  easier  to  determine  experi¬ 
mentally  the  wall  motions  during  instability.  Consequentially,  it  is 
relatively  easy  to  obtain  good  results  in  this  case.  A  typical  set  of 
data,  obtained  from  a  test1  on  an  orthotropic  cylinder,  is  given  in 
Figure  6. 

The  literature  contains  much  material  with  reference  to  the  instability  of 
spheres  and  spherical  caps  under  the  action  of  normal  point  loads,  uniform 
external  pressure  and  combinations  thereof.  Hi  is  has  been  examined  in 
Reference  18. 

19  20 
The  experimental  work  of  Ashwell  and  Evan-Ivanowski,  Cheng,  and  Loo  on 

spherical  caps  subjected  to  a  concentrated  load  at  the  apex  provides 

several  cases  in  which  the  load-displacement  curves  can  be  successfully 

analyzed  by  the  Southwell  method.  The  load-displacement  curve  and  the 

corresponding  Southwell  line  are  given  in  Figure  7  for  Ashwell 's  test  on  a 
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FIGURE  4.  Lund quist -Type  Plot  for  Thin-Walled 

Circular  Cylinder  under  Axial  Compression 
Observation  Point  No  1. 
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FIGURE  5.  Lundquist-Type  Plot  for  Thin-Wailed 

Circular  Cylinder  Under  Axial  Compression 
Observation  Point  No  2. 
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FIGURE  6.  Torsion  Versus  Inward  Wall  Motion  and 
Corresponding  Southwell  Plot  for  Ortho¬ 
tropic  Cylinder. 
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FIGURE  7*  Load -Deflection  Relationship  and  Corres¬ 
ponding  Southwell  Plot  for  Sphere  with 
Point  Load  at  Apex  (Ref  19)  . 
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specimen  with  a  X  of  6.4.  The  corresponding  case  f  >r  the  Fvan-Ivanowski, 
Cheng,  and  Loo  research  (\-6.3)  is  given  in  Figure  6.  The  critical  values 
of  P1  R/Et'3  derived  from  the  plots  are  in  very  close  agreement,  being  3*09 
and  3*1  units,  respectively. 

Although  the  majority  of  the  data  cannot  be  analyzed  in  this  fashion,  it  is 
interesting  to  note  that  those  results  which  can  be  so  treated  appear  to  be 
consistent  and  to  agree  reasonably  well  with  Ashwell's  computed  values.  It 
would  seem,  however,  that  before  a  positive  statement  is  made  in  this  regard, 
another  series  of  tests  should  be  made  with  the  method  of  correlating  via 
load  displacement  firmly  in  mind  from  the  onset.  The  behavior  of  a  spherical 
cap  under  uniform  external  pressurt^has  been  considred  analytically  and 
experimentally  by  Kaplan  and  Fung."'  The  majority  of  their  data  is  amenable 
to  this  method  of  analysis,  and  a  typical  5/p  versa.  6  plot  is  given  in 
Figure  9*  In  Figure  10,  the  deduced  critical  pressures  are  plotted  as  a 
function  of  the  geometric  parameter  \,  which  describes  the  caps. 

It  is  seen  that  all  the  points  are  on  a  smooth  curve.  The  lower  portion  of 
this  curve  corresponds  very  well  with  that  computed  by  Kaplan  and  Fung  as 
an  approximate  solution  to  their  equations. 

The  spherical  cap  under  combined  loading,  a  normal  force  at  the  apex  and  a 
distributec^pressure,  has  been  examined  experimentally  by  Loo  and  Evan- 
Ivanowski.  Many  of  their  load-displacement  curves  can  be  analyzed  by  the 
Southwell  method,  but  there  is  no  theory  available  for  comparison. 

The  behavior  of  complete  spheres  under  external  pressure  loading  is  another 
area  ir.  which  problems  of  correlating  test  data  currently  exist. 

An  experimental  study  was  made  on  a  thin-walled  nickel  sphere  to  check  the 
possibility  of  relating  radial  displacement,  pressure,  and  theoretical 
critical  external  pressure.  The  experiment  is  reported  in  fullgin  Refer¬ 
ence  18.  The  basic  material  had  a  Young's  modulus  of  29*3  X  10  lb/in.  as 
determined  by  routine  material  tests.  This  value  was  confirmed  for 
several  points  on  the  sphere  by  measurement  of  decrease  in  radius  as  a 
function  of  pressure.  The  shell  was  elastically  buckled;*  when  the  point 
of  buckle  inception  had  been  thus  established,  a  test  was  made  using  a 
sensitive  noncontacting  displacement  probe  located  at,  or  very  near,  the 
buckle  center.  The  variation  of  displacement  with  pressure  is  clearly 
seen  in  Figure  11.  There  are  two  properties  of  this  curve  which  are  both 
apparent  and  important.  First,  the  initial  deflection  is  linear.  Second, 
the  final  displacement  curve  is  hyperbolic.  The  hyperbolic  curve  gives 
an  excellent  Southwell  plot  (Figure  12 ),  whose  slope  corresponds  closely 


*The  sphere  was  restrained  again*=l  excessive  motion  by  the  use  of  an 
internal  mandrel.  As  a  result,  on  the  second  and  subsequent  buckling 
tests,  the  load  levels  for  instability  reached  the  same  value  as  on  the 
initial  test. 
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FIGURE  8.  Southwell  Plot  for  Sphere  (Data  from 
Ref  20). 
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FIGURE  9.  Southwell  Plot  for  Spherical  Cap  under 
External  Pressure  (Ref  2l) . 


GEOMETRIC  PARAMETER  A 


FIGURE  10.  Critical  Pressure  Versus  \  as  Deter¬ 
mined  by  Southwell  Plots  (ref  2l). 
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FIGURE  11.  Pressure-Displacement  Relation  and 

Corresponding  Southwell  Plot  for  Nickel 
Sphere  under  External  Pressure. 


FIGURE  12.  Southwell  Plot  for  Sphere  External 
Pressure  (Data  from  Ref  23). 
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to  the  critical  load  computed  on  classic  theory  when  he  modulus  has  the 
value  indicated  by  the  initial  slope. 

The  data  presented  by  Thornpson‘-'>have  also  been  examined  in  the  same  manner 
The  conclusions  reached  were  identical.  The  Southwell  plot  predicts  the 
classic  critical  load  for  this  sphere,  also.  It  should  he  noted  that  in 
this  case,  too,  the  buckle  generating  defect  was  located  at  a  point  which 
was  characterized  by  an  apparent  reduced  modulus.  The  modulus  value  used 
in  the  data  correlation  was  that  which  corresponded  to  the  initial  linear 
portion  of  the  load-displacement  curve. 
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CONCLUSIONS 


The  experimental  and  theoretical  study  presented  in  this  report  concludes 
a  broad  analysis  which  has  been  conducted  on  the  range  of  applicability  of 
the  Southwell  plot  as  a  means  of  data  interpretation  for  tests  on  shells. 

The  results  presented,  when  taken  in  conjunction  with  those  given  in  a 
previous  report  for  plate  structures,  provide  evidence  that  the  method  is  of 
absolute  generality.  However,  it  may  be  reasoned  in  the  case  of  structures 
in  which  many  buckles  can  be  formed,  that  the  true  process  of  buckuing  may 
be  controlled  by  the  individual  character  of  the  structures.  In  this  case, 
imperfection  could  influence  local  boundary  values  and  would  thus  show  as 
an  apparent  lowering  or  an  elevating  in  critical  load  as  given  by  the  South- 
well  plot.  If  this  i.  so,  we  should  anticipate  that  the  distribution  of 
Southwell  values  will  be  essentially  normal  if  many  observations  are  made  at 
devious  points  over  the  body. 

With  this  process,  then,  we  feel  confident  in  asserting  that  theories  for 
perfect  todies  which  give  estimates  not  in  agreement  with  experimental  values 
derived  on  the  basis  of  Southwell- type  plots  must  be  considered  inadequate, 
incomplete,  or  in  error. 
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Thie  report  demonstrates,  from  experimental  studies,  that  the  classic  critical  } 
loads  computed  from  small  displacement  theory  are  correct  for  perfect  shells  under  j 
ideal  loading  conditions.  Hie  method  used  is  an  implicit  rather  than  an  explicit  * 
one.  It  makes  use  of  the  fact  that  for  a  realistic  structure,  the  elastic 
deformation  can  be  associated  with  the  initial  displacement  from  ideal  form,  the 
load  which  produces  the  motion,  and  the  classic  load  for  the  structure  by  the 
hyperbolic  expression  k 
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By  choosing  the  variables  to  be  6/P  and  6/  the  relationship  can  be  presented  in  the 
form  of  a  straight  line  whose  slope  corresponds  to  the  critical  load  for  the  ideal 
case.  This  is,  in  essence,  the  method  developed  by  Arytor  and  Perry,  in  1869,  to 
analyze  column  data.  The  generalization  to  include  other  structures  was  foreseen 
by  Southwell  in  his  classic  paper  of  1932,  although  he  offered  no  proof.  In  effect 
a  general  proof  exists  in  the  theory  of  elastic  stability  as  presented  by  Westtr- 
gaard  in  1922.  However,  the  practical  application  to  shells  has,  until  now,  not 
been  made. 

By  using  this  technique,  the  behavior  of  cylinders  under  the  action  of  external 
pressure,  torsion,  and  axial  load  and  of  the  instability  of  spheres  and  spherical 
caps  under  the  action  of  external  pressure  is  examined.  In  the  main,  the  analyses 
are  conducted  on  experimental  data  already  published.  However,  there  are  notable 
exceptions,  the  cylindrical  shell  under  axial  load  a  case  In  point. 
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